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We study a random neighbor version of the Bak-Sneppen model, where "nearest neighbors" are 
chosen according to a probabihty distribution decaying as a power-law of the distance from the 
active site, P{x) ~ |x — Xacl'"^ ■ AH the exponents characterizing the self-organized critical state of 
this model depend on the exponent oj. As a; —> 1 we recover the usual random nearest neighbor 
version of the model. The pattern of results obtained for a range of values of ui is also compatible 
with the results of simulations of the original BS model in high dimensions. Moreover, our results 
suggest a critical dimension dc = 6 for the Bak-Sneppen model, in contrast with previous claims. 
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Since its introduction, the Bak-Sneppen (BS) model ||] 
has had much success as perhaps the simplest and yet 
non-trivial self-organized critical (SOC) extremal model. 
Understanding its behavior is therefore very important 
to get an insight in the behavior of other SOC extremal 
models §. 

The BS model is easily defined: To each site i on a 
hypercubic lattice in d dimensions is assigned a random 
variable fi taken from a probability distribution p{f), 
say, uniform in [0,1]. Then at each time-step the site 
i with the smallest fi is chosen (it is called the active 
site) , and its variable and the variables of its 2d nearest- 
neighbors are updated taking them from p{f). As a re- 
sult of this dynamics, the system organizes in a station- 
ary state where almost all the variables fi are above a 
threshold fc- Moreover, in this state, the dynamics of 
the model has self-similar features: each update of the 
minimum variable triggers a local avalanche of updates; 
the time durations of the avalanches obey a power-law 
distributions characterized by an exponent r. The num- 
ber of sites touched by an avalanche of duration t grows 
like f^. Also the first return times (defined as the times 
between two successive returns of the activity to the same 
site) and the all return times (the times between the first 
passage of the activity on a site and any successive re- 
turn to the same site) are power-law distributed, with 
exponents r/ and Ta respectively. 

Not all of the above exponents are independent. In- 
deed it is possible to show, from renewal theory, that 
Ta + T/ = 2 if Tq < 1, Ta = T/ if Ta > 1 ||] . Recently, 
a non trivial relation between t and /x has been unveiled 
in MM exploiting the hierarchical structure of the up- 
date avalanches (each avalanche is made up of smaller 
avalanches, and so on down to the microscopic scale). 
Both relations are satisfied for d = 1 with t ~ 1.07 and 
^l ~ 0.42, Ta ^ 0.42 and t/ - 1.58 [^. 

The only known exactly solved version of the BS model 
is the Random Nearest-Neighbor (RNN) model: there 



"nearest neighbors" are chosen at random over the lat- 
tice Q . As a result geometric correlations typical of low 
dimensions are lost, and the RNN can be considered as 
a Mean Field version of the BS model. The exponents of 
the RNN model are known to be t = Tq = t/ = 3/2 
and /i = 1 ra. In particular these exponents satisfy 
both Ta = Tf and the relation between t and /i proposed 
in Q. In js) this relation has been carefully studied, 
and it has been "graphically" explicited (see Fig.|l|). The 
knowledge of the two exponent relations has given con- 
fidence in high-dimensional simulations S whose main 
result is that the upper critical dimension of the model 
is due = 8 (the upper critical dimension is the dimension 
above which the exponents should take the RNN values) . 
This conclusion is at odds with previous claims, based on 
analogies of the BS model with directed percolation, that 
due should be 4. A further result from S is the presence 
of two regimes as the dimension d of the system increases: 
for d < 3 the model is recurrent (tq < 1; in random walk 
theory recurrence means that every site of the lattice 
is touched an infinite number of times with probability 
one), whereas for d > 4 the model is transient (tq > 1; 
transience means that there is a finite probability, smaller 
than 1, that a site will be touched by the process), yet 
non trivial (i.e., different from the RNN model) as long 
as d < 8. Actually, before g], a further exponent rela- 
tion was believed to hold, namely ji — Ta- In [p[ this 
relation is shown not to hold for d > 2. Indeed, since 
/i < 1 always and Ta{RNN) — 3/2, it is straightforward 
to conclude that at least as soon as Tq > 1, Ta ^ /i. The 
dimensionality d = 3, with Ta < 1 and ^l ^ Ta, can be 
therefore considered as representative of a further regime 
within the recurrent one. 

The BS model shows therefore an extremely rich be- 
havior changing the dimensionality d of the system. Yet, 
high-dimensional simulations are always susceptible of 
strong finite-size corrections, and the good convergence 
of the results is difficult to prove. In this Communication 



we propose a new way to interpolate between the d — 1 
and the RNN models: The "nearest neighbors" of the 
active site Xac are chosen at random over the lattice, but 
with a probability that is a power-law decreasing function 
of the distance from it 
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As we will show, varying lo we find the same behavioral 
pattern as found in |q| varying d. We name this model 
the Levy-Random-Nearest-Neighbor model (LRNN; here 
the use of the word Levy is somehow an abuse since we 
use also w > 3). 

As a loose analogy, we recall that the same idea has 
been applied also to the d = 1 Ising model with interac- 
tions decaying as (0), and indeed it has been found that, 
varying w it is possible to go from the d — 1 model to 
mean- field like results Q. 

Simulations are performed over Id lattices of up to 2^^ 
sites, with growing sizes showing stability of the expo- 
nents. 

In Figjl] we show the r avalanche exponents for differ- 
ent values of to plotted against the corresponding values 
of jj. All the fx/r pairs nicely satisfy the exponent rela- 
tion between the two exponents obtained in ||]. This is 
a first important check of the consistency of our simula- 
tions and of the exponent relation. 

In FigJa we plot the Tq, t/ and /i exponents for dif- 
ferent values of to from ui — 3 to ui = 1^ (this extreme 
value is not shown since simulations become extremely 
difficult due to the non normalizability of the distribu- 
tion if u! — 1). Many important aspects of the model can 
be discussed looking at Fig.0. We find that the exponent 
relation between Ta and t/ is satisfied both when Ta < I 
and when Ta > I. Therefore this result and Fig.0 confirm 
the validity of the two already known exponent relations. 
Moreover, we see that the exponents tend to their RNN 
values as u! —^ 1. This result should have been expected. 
Indeed the probability distribution (|l|) is normalizable in 
the thermodynamic limit only as long as w > 1; when 
Lu < 1 then the normalization is ruled by the length of 
the lattice 
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that diverges in the limit of infinite lattice size L. There- 
fore the distribution P{x), properly normalized, degen- 
erates to 0, just as in the RNN case, where the normal- 
ization is \/L (case oj = 0). 

The opposite case, u) — > oo, is also intriguing. Indeed 
we could naively expect the d — 1 limit to be recovered 
when a; > 3. In that case the average distance of the 
"neighbors" and its variance are finite, just as for near- 
est neighbors. We would thus expect the w > 3 case 
to belong to the same universality class as the original 
BS model, but this conclusion is not correct. In order 



to shed light on this problem, we also performed simula- 
tions taking the neighbors according to a distribution ex- 
ponentially decreasing with respect to the distance from 
the active site Xac- In this case, instead, we nicely recover 
the known d = 1 exponents of the BS model. We con- 
clude therefore that the presence of diverging moments of 
order higher than two drives the system out of its nearest 
neighbor fixed point. The latter holds instead whenever 
the distribution of the random neighbors has all its mo- 
ments finite. This result, although non trivial, is not new 
in (annealed or quenched) disordered systems. For exam- 
ple, it has been shown that diverging moments of order 
higher than two in the disorder distribution can change 
the universality class of directed polymers in random en- 
vironments, and of the related Kardar-Parisi-Zhang sur- 
face growth equation ||l^ . 

It is well known in classical random walk theory that 
random walks with a jump probability distribution with 
finite variance belong to the Gaussian universality class. 
Random walks with infinite higher moments have a mi- 
croscopic structure that is different from a Gaussian one, 
mainly made of clusters of points seldom separated by 
long jumps. On larger and larger length scales, this clus- 
ter structure disappears, and the walks "renormalize" to 
Gaussian ones. Yet, a main difference between a simple 
random walk and the BS model is the presence of mem- 
ory effects. Indeed, any time a site is chosen (either as 
the active one or as one of its neighbors), any memory 
of its previous updates is lost. Therefore, the interac- 
tion between the geometric structure given by the choice 
of the neighbors according to P{x) and the updates of 
the corresponding variables can give rise to non-trivial 
effects. 

As already mentioned above, the geometrical fractal 
dimension Df of the avalanches is another quantity of 
interest. It is possible to relate Df to the other expo- 
nents of the model. Indeed, we recall that the fractal 
dimension relates the volume N of the avalanche (that 
is, the number of sites touched by the avalanche) with 
the typical size R of the avalanche, as iV ~ R^f . On the 
other hand, N scales with the duration of the avalanche 
a.s N ^ f^. The relation between the typical size of the 
avalanche R, and its duration t is given by i ~ i?^, z be- 
ing the dynamical exponent of the model. Then we find 
Df — zfi. 

In principle, in order to determine z it is possible to use 
the all return time distribution Pa{t). At long times on 
a d— dimensional lattice oi N = L'^ sites, Pa{t) flattens. 
W e can write a scaling form for Pa [t) as 
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with f{x) ^ const when x —>■ and f{x) ~ a;^" when 
a; — > oo. In this second case we find that Pa{t, L) = L"'* 
(roughly speaking, as soon as every site of the lattice has 



been touched by an avalanche, they have all the same 

probability to be chosen), from which we have z = dJTa- 

Then we can write an expression for the fractal dimen- 
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This expression holds for the high dimensional simu- 
lations of [pi: when avalanches are compact objects 
{d = 1,2), fj, = Ta- Then, fi j^ Ta and avalanches 
become fractal objects. Indeed (0) approximates well 
the data given in g for d = 3 (va = 0.92, fi = 0.85, 
Df = 2.6 ~ 3^/Ta = 2.77...) and for rf = 4 {Ta = 1.15, 
/i = 0.92, Df = 3.3 ~ Afi/Ta == 3.2). As a byprod- 
uct, we find that (j^ suggests an upper critical dimen- 
sion for the Bak-Sneppen model due = 6: indeed with 
/^ = 1 and r = 3/2 (the "mean-field", RNN, values of 
the exponents) we find Df{d = 6) = 4, that is indeed the 
predicted avalanche fractal dimension in the RNN limit. 
due = 6 is at odds with what stated in H , where due = 8 
was suggested by numerical simulations, but also with 
Pi where due = 4 was claimed based on analogies with 
directed percolation. 

Although (Q) seems to hold in the high dimensional 
case, it does not fit the numerical results in the present 
LRNN approach. The reason is that, whereas in the high 
dimensional case there is a single relation between time 
and space, namely t ~ L''/'^", in the LRNN case there is 
a further relation, the usual Levy random walk law t ~ 
R^^^ . We measure the fractal dimension of avalanches 
using the distance between the rightmost and leftmost 
touched sites as a measure of i?, and this corresponds 
to z ^ Lo — \. Indeed, as it can be seen from FigJ3, 
Df = /i(u; — 1) approximates very well the measured 
fractal dimensions. 

As noted above about the fractal dimension of high- 
dimensional avalanches, we observed that ^ = Ta corre- 
sponds to compact avalanches. In the LRNN case we see 
from Fig.|| that indeed /i 7^ Tq for cj < 2, even if the 
fractal dimension Df < 1 already for w < 3. We can try 
to understand this result remembering that D f is related 
to the random walk exponent z = uj ~ 1: although z is 
different from its Gaussian value z = 2 as soon as w < 3, 
a Levy random walk with w > 2 is still compact, and 
so is the structure built by a choice of neighbors accord- 
ing to (nl). Only when lo < 2 such a structure becomes 
genuinely fractal, and [i ^ Ta- 

The possibility to obtain the exponent D f from /i and 
Ta in high dimensions, and from /i and 10 in the LRNN 
version, suggests that indeed there are at most two inde- 
pendent exponents in the model, namely /i and Tq. The 
LRNN model suggests the presence of a further (although 
non trivial) exponent relation. Indeed, in Fig.0 we show 
the values of the r^ exponent as a function of the cor- 
responding /i exponent for different values of w and for 
different dimensions. As it can be seen, the agreement is 



good, suggesting that the knowledge of /x (or of Tq) is suf- 
ficient to know all the other exponents through relations 
that, as for the [iJt one, could be highly non-trivial. 

In conclusion, we have introduced a modification of the 
Bak-Sneppen model where the neighbors of the active site 
are chosen at random over the lattice with a probability 
that decreases like a power-law of the distance from the 
active site, with an exponent lo. As a result we find that 
the characteristic exponents of the model interpolate be- 
tween the d—\ limit (u) = c») and the mean field (RNN) 
limit (w < 1). In particular, we verify that the known 
exponent relations hold for this model too. Moreover, we 
find and verify an exponent relation for the fractal dimen- 
sion of the avalanches, Df — (w — l)/i. As a byproduct we 
obtain a relation between Df and /i and Ta also in high 
dimensions, fitting well the present numerical results up 
to d = 4 g] , although it suggests an upper critical dimen- 
sion due — 6 (and not due = 4 or due = 8 as previously 
believed). More accurate numerical simulations in high 
dimensions are therefore needed. The relevance of the re- 
sults reported in this Communication is manifold: they 
can be looked at as an interesting modification of the 
Bak-Sneppen model, but their full importance emerges 
when compared to the high dimensional results of [gj. 
Indeed, they lead us to propose a new value of the upper 
critical dimension of the model, namely due = 6, and to 
conjecture the existence of a still undiscovered exponent 
relation between fj, and Ta, reducing therefore the num- 
ber of independent exponents to one in any dimension. 
This results show therefore that there is still some way 
to go before a full and satisfying understanding of the BS 
model is achieved. 

The authors thank F. Slanina for useful discussions. R. 
Cafiero and P. De Los Rios aknowledge financial support 
under the European network project FMRXCT980183. 



[1] P. Bak and K. Sneppen, Phys. Rev. Lett. 71, 4083 (1993). 

[2] M. Paczuski, S. Maslov, and P. Bak, Phys. Rev. E 53, 
414 (1996). 

[3] M. Fisher, in Fundamental Problems in Statistical Me- 
chanics VI, ed. E.G.D. Cohen, North-Holland (Amster- 
dam 1984). 

[4] S. Maslov, Phys. Rev. Lett. 77, 1182 (1996). 

[5] M. Marsili, P. De Los Rios and S. Maslov, Phys. Rev. 
Lett. 80, 1457 (1998). 

[6] H. Flyvbjerg, K. Sneppen, and P. Bak, Phys. Rev. Lett. 
45, 4087 (1993). 

[7] J. de Boer, B. Derrida, H. Flyvbjerg, A.D. Jackson, and 
T. Wettig, Phys. Rev. Lett. 73, 906 (1994). 

[8] P. De Los Rios, M. Marsili and M. Vendruscolo, Phys. 
Rev. Lett. 80, 5746 (1998). 

[9] S.A. Cannas, Phys. Rev. B 52, 3034 (1995). 



[10] T. Halpin-Healy and Y.-C. Zhang, Phys. Reports 254, 
216 (1995). 



1.60 
1.50 



D, 



0.6 



0.4 - 



0.2 



1 




1.5 



2 
CO 



2.5 



FIG. 3. Avalanche fractal dimension D[ for different values 
of Lo (circles) and as from Df = (cj — l)p (squares). 



FIG. 1. Avalanche exponent r vs. // for different Levy 
exponents w. The dashed line is the exact relation between 
the two exponents as from [5] . The data from the simulations 
and the exact relation are in excellent agreement. Error bars 
are comparable with the symbol dimensions. 
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FIG. 2. Values of the different relevant exponents, /x, Ta 
and T[, as a function of the Levy exponent w. 
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FIG. 4. All return time exponent Ta vs. /i for different 
values of u) (squares) and for different dimensions, from d = 1 
to d = 8 (circles). 



